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t>- ■ Abstract 
00 

^ The energy loss effect in nuclear matter is another nuclear effect apart from 

the nuclear effects on the parton distribution as in deep inelastic scattering 
process. The quark energy loss can be measured best by the nuclear 
dependence of the high energy nuclear Drell-Yan process. By means of 
three kinds of quark energy loss parameterizations given in literature and 
the nuclear parton distribution extracted only with lepton-nucleus deep 
inelastic scattering experimental data, measured Drell-Yan production 
cross sections are analyzed for 800GeV proton incident on a variety of 
nuclear targets from FNAL E866. It is shown that our results with 
considering the energy loss effect are much different from these of the 
FNAL E866 who analysis the experimental data with the nuclear parton 
distribution functions obtained by using the deep inelastic 1A collisions 
and pA nuclear Drell-Yan data . Considering the existence of energy loss 
effect in Drell-Yan lepton pairs production, we suggest that the extraction 
of nuclear parton distribution functions should not include Drell-Yan 
experimental data. 
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Understanding the initial stages of ultrarelativistic heavy ion collisions is of 
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utmost importance in order to understand the outcome of the high energy heavy 
ion experiments, such as the BNL relativistic heavy ion collider(RHIC) and CERN 
large hadron collider (LHC). Understanding the modifications of the parton distri- 
bution functions and the parton energy loss in nuclei should be the first important 
step towards pinning down the initial conditions of a heavy-ion collision and un- 
derstanding of J/ip production which is required if it is to be used as a signal for 
the quark-gluon plasma in relativistic heavy ion collisions. 

The high energy inelastic hadron-nucleus collisions has been studied for many 
decades by both the nuclear and particle physics communities M . By means of the 
nuclei, we can study the space-time development of the strong interaction during 
its early stages, which is inaccessible between individual hadrons. In high energy 
inelastic hadron-nucleus scattering, the projectile rarely retains a major fraction 
of its momentum after traversing the nucleus. Rather, its momentum is shared by 
several produced particles, which form a hadron jet in the forward direction. The 
classical description of this phenomena is that the projectile suffers multiple colli- 
sions and repeated energy loss in the nuclear matter. In other words,each quark or 
gluon in the projectile can loss a finite fraction of its energy in the nuclear target 
due to QCD bremsstrahlungt 2 ' . The Drell-Yan reaction^ on nuclear targets pro- 
vides, in particular, the possibility of probing the propagation of quark through 
nuclear matter, with the produced lepton pair carrying away the desired informa- 
tion on the projectile quark after it has travelled in the nucleus. Only initial-state 
interactions are important in Drell-Yan process since the dimuon in the final state 
does not interact strongly with the partons in the nuclei. This makes Drell-Yan 
scattering an ideal tool to study energy loss' 4 l In addition, Drell-Yan reaction is 
closely related to DIS of leptons , but unlike DIS ,it is directly sensitive to anti- 
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quark contributions in target parton distributions. When DIS on nuclei occurs at 
x < 0.08, where x is the parton momentum fraction, the cross section per nucleon 
decreases with increasing nucleon number A due to shadowing' 5 ' . Shadowing should 
also occur in Drell-Yan dimuon production at small X2, the momentum fraction of 
the target parton, and theoretical calculations indicate that shadowing in the DIS 
and Drell-Yan reactions has a common origin '^.Therefore, shadowing and initial 
state partonic energy loss are processes that occur in the proton-induced Drell-Yan 
reaction on nuclei. 

In order to describe the modification of the initial state parton distributions 
in nucleus , a variety of approach to this question exist in the literature^. Re- 
centlythere were two trials without model dependence to obtain nuclear parton 
distributions from the existing world experimental data. In 1999, Eskola, Kolhi- 
nen, Ruuskanen and Salgado(EKRS)® suggested a set of nuclear parton distri- 
butions ,which are studied within a framework of the DGLAP evolution . The 
measurements of FA/F® in deep inelastic I A collisions, and Drell-Yan dilepton 
cross sections measured in pA collisions were used as constraints. The kinematical 
ranges are 10~ 6 < x < 1 and 2.25GeV 2 < Q 2 < 10 A GeV 2 for nuclei from deuteron 
to heavy ones .With the nuclear parton distributions , the calculated results agreed 
very well with the relative EMC and Fermilab E772 experimental data' 9 ' . In 2001, 
Hirai,Komano and Miyama(HKM)[ 10 l proposed two types of nuclear parton distri- 
butions which were obtained by quadratic and cubic type analysis, and determined 
by a x 2 global analysis of existing experimental data on nuclear structure functions 
without including the proton-nucleus Drell-Yan process .The kinematical ranges 
covered 10~ 9 < x < 1 and lGeV 2 < Q 2 < 10 5 GeV 2 for nuclei from deuteron to 
heavy ones . As a result, they obtained reasonable fit to the measured experimental 
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data of F 2 .In this work , we will use both parameterizations and investigate the 
nuclear dependence of the Drell-Yan process. 

In order to research the partonic energy loss in nuclei, Fermilab Experi- 
ment866(E866) ^ performed the precise measurement of the ratios of the Drell- 
Yan cross section per nucleon for an 800GeV proton beam incident on Be, Fe and 
W target at larger values of x±,the momentum fraction of the beam parton, larger 
values of if(~ %i — #2), and smaller values of X2 than reached by the previous ex- 
periment, Fermilab E772' 9 '. The extended kinematic coverage of E866 significantly 
increases its sensitivity to energy loss and shadowing. This is the first experiment 
on the energy loss of quark passing through a cold nucleus. 

After the E866 experimental data was reported,several groups have given 
their theoretical analysis of the data^ 12,13,14 ]. In previous report^ 14 ' , by means of 
EKRS and HKM nuclear parton distribution functions, we investigated the Drell- 
Yan production cross section ratios for 800GeV proton incident on a variety of nu- 
clear targets in the framework of Glauber model. We found that the theoretical re- 
sults with energy loss are in good agreement with the Fermilab E866 experiment by 
means of HKM nuclear parton distributions. However, the calculated results with- 
out energy loss can give good fits by using EKRS nuclear parton distribution func- 
tions. In this letter, the nuclear dependence of the pA Drell-Yan production cross 
sections are studied by combining the quark energy loss parametrizationl 15,18 ' 19 ' 
given in literature with the EKRS and HKM nuclear parton distribution. 

In the Drell-Yan process,the leading-order contribution is quark-antiquark 
annihilation into a lepton pair. The annihilation cross section can be obtained 
from the e + e~ — > cross section by including the color factor | with the 
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charge e'j for the quark of flavor /. 

da 8na 2 ~ 2 

where = (xia^s) 1 / 2 ^ the center of mass system (cm. system) energy of qq col- 
lision ,xi(resp.a;2)is the momentum fraction carried by the projectile (resp. target) 
parton, y/s is the center of mass energy of the hadronic collision, and M is the 
invariant mass of the produced dimuon. The hadronic Drell-Yan differential cross 
section is then obtained from the convolution of the above partonic cross section 
with the quark distributions in the beam and in the target : 

where K is the high-order QCD correction, a is the fine-structure constant, the 
sum is carried out over the light flavor / = u,d,s, and q [ f A \x) and <ff A \x) are 
the quark and anti-quark distributions in the proton (nucleon in the nucleus A). 
In order to obtain the x\ dependence of Drell-Yan production, we shall deal in the 
following with the single differential cross section, 

^^E^/f^^N+^l^)], (3) 

where the integration over the dimuon mass is performed in the range given from 
E866 experiment. 

Now let us take into account of the energy loss of the fast quarks moving 
through the cold nuclei. There were three expressions for the average change in in- 
cident parton momentum fraction in literature .The first was given by Gavin and 
Milana^ 15 ' in 1992. They analyzed the depletion at high Feynman x in large rela- 
tive to small nuclear targets from measurements of the Drell-Yan process and of 
charmonium and bottomonium production at CERNl 17 ] and Fermilab E772I 9 ' 16 ! . 



They proposed that energy loss due to multiple parton scattering can consistently 
explain the high-Xi? depletion in the Drell-Yan process and J/ip production. The 
following expression was assumed for the average change in incident parton mo- 
mentum fraction : 

Ax 1 = K1X1A 1 / 3 , (4) 

where the factor K\ may have a Q 2 dependence. From comparison to the exper- 
imental data and neglecting shadowing, an upper bound on dE/dz-the incident 
quark energy loss per unit length — was obtained, the value is 1.5GeV/ fm.This 
result was questioned by Brodsky and Hoyer' 18 ! ,who argued that the time scale 
for gluon bremsstrahlung need to be taken into account. By analogy with the pho- 
ton bremsstrahlung process , they derived the quantum mechanical bound on the 
amount of radiative energy loss suffered by high energy quarks and gluons in nuclear 
matter, 

Axi ~ — A 1/3 , (5) 
s 

where s denotes the square of the proton-nucleus center of mass energy. The 
formulation developed by Brodsky and Hoyer was extended by Baier et al.' 19 ' They 
studied the medium-induced p^-broadening and induced gluon radiation spectrum 
of a high energy quark or gluon traversing a large nucleus. Multiple scattering of 
the high energy parton in nuclei was treated in the Glauber approximation. Baier 
et al. found that the partonic energy loss depended on a characteristic length 
and the broadening of the squared transverse momentum of the parton. For finite 
nuclei,both factors vary as A 1 / 3 , so they predicted 

Axi « —A 2/3 . (6) 

s 

Using these energy loss expressions , we can obtain empirical values for the «'s by 
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fitting the FNAL E866 Drell-Yan cross section ratios versus the incident proton's 
momentum fraction and dimuon effective mass. 

If considering the quark energy loss in nuclei, the incident quark momentum 
fraction can be shifted from x[ — x± + Ax± to x\ at the point of fusion. Combining 
the shadowing with initial state energy loss ,the production cross section in pA 
Drell-Yan process can be written as 

In order to pin down quark energy loss by comparing with the experimental data 
from E866 collaboration! 11 ! ,we introduce the nuclear Drell-Yan ratios as: 

/ x da p ~ Al ,da p - M 

= —I— (8) 

The integral range on M is determined according to the E866 experimental kine- 
matic region. In our theoretical analysis ,x 2 is calculated with the Drell-Yan differ- 
ential cross section rations Ra 1 /a 2 as 

/ r>data T>theo \2 

2 = V- [^Ai/Mj n A 1 /A 2 ,j) (Q s 
A 2-*i ( perr \2 ' v ' 

3 y U A 1 /A 2 ,j) 

where the experimental error is given by systematic errors as R e A[/A 2 ,ji an< ^ ^-M^ji 
R t A^/A 2 j) indicates the experimental data (theoretical values ) for the ratioi?Ai/A 2 - 
Taking advantage of the EKRS' 8 ' nuclear parton distribution functions and 
Ki = 0,k 2 = 0,/t 3 = (without energy loss effects), the obtained x 2 value is 
X 2 = 51 .4 for the 56 total data points. The x 2 P er degrees of freedom is given 
by x 2 /d.o.f. = 0.918. It is apparent that theoretical results without energy loss 
effects agree very well with the E866 experimental data , which results from EKRS 
parametrization of nuclear parton distributions studied with including the E772 
Drell-Yan data ^ 16 1 . Because of employing the EKRS global phenomenological 
analysis, the E866 analysis missed the energy loss effect of quarks in nuclei. 



In what follows, we consider also combining HKM cubic type of nuclear parton 
distribution' 10 ^ with the three quark energy loss parameterizations. With n\ = 
0, K2 = 0, K3 = (without energy loss effects), the obtained x 2 P er degrees of freedom 
is x 2 /d.o.f. = 2.516. With m = 0.0036, k 2 = 3.45, k 3 = 0.46(with energy loss 
effects), the obtained x 2 P er degrees of freedom are x 2 /d.o.f. = 1.00048, 1.022, 1.053, 
respectively. The results given by HKM quadratic type are nearly the same as these 
above. As an example, the calculated results with energy loss expression Eq.(5)are 
shown in Fig.l and Fig. 2. which is the Drell-Yan cross section ratios for Fe to 
Be and W to Be as functions of Xi for various interval of M,respectively.The solid 
curves are the ratios with only the nuclear effect on the parton distribution as in 
DIS scattering process, and the dotted curves correspond to an energy loss effect: 
k 2 = 3A5GeV 2 with nuclear effect on structure function .From comparison with the 
experimental data , it is found that our theoretical results with energy loss effect are 
in good agreement with the Fermilab E866. However, the values of Ki, k 2) k 3 are very 
different from those by the E866 Collaboration' 11 !, which results from employing 
EKRS nuclear parton distributions. The EKRS analysis itself included the Drell- 
Yan data from E772,with the presumption that the low-a^ nuclear dependence 
arose entirely from shadowing. This clearly introduced an inconsistency into the 
E866 search for energy loss. Therefore, we can conclude that the extraction of 
nuclear parton distribution functions should not include Drell-Yan experimental 
data because of the energy loss of the incoming partons. 

In Eq.(4), the average loss of projectile parton vary with the variable x±. 
It is difficult to obtain a constant energy loss dE / dz(GeV / fm).K 3 = 0.46 shows 
us that the observed energy loss of the incident quarks is AE < 0.21GeV/ fm 2 x 
L 2 .k 2 = 3.45 in Eq.(5)indicates that the incident quarks lose energy with a con- 
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stant rate dE/dz = 2.0AGeV/fm if < L >= 3/4(1. 2A 1 / 3 )fm is used as the average 
path length. The value of dE/dz is not obviously consistent with that given by 
M.B.Johnson et aU 12 l, which is -dE/dx = 3.14 ± 0.53GeV//m.M.B. Johnson et al. 
examined the effect of initial state energy loss on the Drell-Yan cross section ratios 
from E866 by employing a new formulation of the Drell-Yan process in the rest 
frame of nucleus.lt is obvious that an intuitive explanation of the quark energy 
loss depends on the reference frame of the target nucleus. It is well-known that 
the physical observables can not depend on our choice of reference frame. There- 
fore,further theoretical and experimental studies are needed to determine the value 
of incident quark energy loss per unit length, and verify whether the energy loss is 
linear or quadratic with the path length. 

With considering the existence of energy loss effect in Drell-Yan lepton pairs 
production, we suggest that the extraction of nuclear parton distribution functions 
should not employ pA Drell-Yan experimental data. Although there are currently 
abundant data on electron and moun deep inelastic scattering, it is difficult to 
determine valence quark distributions in the small x region and the anti-quark 
distributions .Only valence quark distributions in medium x region can be relatively 
well determined. It is well considered that the precise nuclear parton distributions 
must be known in order to calculate cross sections of high energy nuclear reactions 
accurately and find a signature of quark-gluon plasma in high energy heavy-ion 
reactions. We suggest using precise neutrino scattering experimental data ,which 
can provide a good method for measuring the F 2 (x, Q 2 ) and xF 3 (x,Q 2 ) structure 
functions. Using the average of xFg A (x, Q 2 )and xFg A (x, Q 2 ) , the valence quark 
distribution functions can be well pinned down' 20 '. Combining the lepton inelastic 
scattering data with the neutrino scattering experiments, valence quark and anti- 
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quark distribution functions will be obtained in the future, which makes us good 
understanding the energy loss effect in high energy nuclear collisions. 

In summary, we have made a leading-order analysis of E866 data in nuclei by 
taking into account of the energy loss effect of fast quarks . Our theoretical results 
with quark energy loss are in good agreement with the Fermilab E866 experiment 
by means of the parametrization of nuclear parton distributions studied without 
nuclear Drell-Yan process, which is the same as that in our previous work^.we 
suggest that the extraction of nuclear parton distribution functions should not em- 
ploy pA Drell-Yan experimental data. In order to determine the amount of energy 
loss and the appropriate path length,we desire to operate precise measurements 
of the nuclear dependent of Drell-Yan production at the Fermilab Main Injec- 
tor (FMI,120GeV proton beam)! 21 ! and the Japan Hadron Facility (JHF, 50 GeV pro- 
ton beam) ^ , where shadowing effect disappears and energy loss effect of fast quarks 
could provide the dominant nuclear dependence. 
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Figure caption 

Fig.l The nuclear Drell-Yan cross section ratios Ra 1 /a 2 { x i) on F e to Be for various 
intervals M. Solid curves correspond to nuclear effect on structure function . Dotted 
curves show the combination of shadowing and energy loss effect with HKM cubic 
type of nuclear parton distributions .The experimental data are taken from the 
E866[ll]. 

Fig. 2 The nuclear Drell-Yan cross section ratios Ra 1 /a 2 { x i) on W to Be for 
various intervals M. The comments is the same as Fig.l 
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